


Institutional Archive of the Naval Postgraduate School 





Calhoun: The NPS Institutional Archive 
DSpace Repository 


Theses and Dissertations l. Thesis and Dissertation Collection, all items 


1991-03 


Structural considerations for aircraft payload 
modification-P-3H zero fuel weight increase. 


Culpepper, Steven Drew 


Monterey, California. Naval Postgraduate School 


http://hdl.handle.net/10945/28567 


Downloaded from NPS Archive: Calhoun 


Calhoun is the Naval Postgraduate School's public access digital repository for 


| К D U DLEY research materials and institutional publications created by the NPS community. 
| Calhoun is named for Professor of Mathematics Guy K. Calhoun, NPS's first 
th 
KNOX appointed — and published — scholarly author. 


WU LIBRARY Dudley Knox Library / Naval Postgraduate School 


411 Dyer Road / 1 University Circle 
Monterey, California USA 93943 








http://www.nps.edu/library 



















































































































































































































































































































































































































































eh 
by = Gu. 
` Pree ét 
PT trn 
ы Mur ch ab e Prem 
PIAL 0 26%. 3 21979 ЫЛА тр 4 m 
met ry "T? AS 5.4. løðir) Аларда 
DINE DIC UT IER Н ГЕ fora r TA v T. 24 
er tre kar d AS PUTT OEY Ин T аб. арт 
AY Dy eo ул P rg арно ^ff ZZ = 
AILE F T ҚЫ, AÐ у а reden) ан HT OUT EL Eee 
ИИ рт. D ETATS d | 
ا‎ N KE Eon ا‎ R C e ENE E EE e 577 
L | шн ФА Ру 204. 83. р str pus % нн M ы 
pr ai apr M HE r „ә М ru B D 4 A ЕУ THEE Tr: en же Лы и элй 
u 4 ы” EYT. EN gh n e. oe TT 177 P n PTS 
ikk " ne Ty mE EL КЕ HÁ ы Pee v 
i TN uw y VI 1 ЕА oum Ur UIA ee RECTE Fk DE 
я || ^ A i m 154 2. ЕРТ 4 Vir I Py HARI s er M “+ XU s^ КШ 1 prie GEO ШЕ т ا و‎ к 
Ld ç к-г ЕСЕ j 4.14. < و‎ 484 и) m ЖАУ hý hbo ӨРҮ Hd M" уч uv. бе! a 4445 
4 ^e Ae "ro^ Bb. KEEN A t Ap f HAP, La FY и! res На e Ha Жтт Pr A „л b 
ха шә, €. M [epa Pu NA MEE AEE KEE To A feb ile УИ x RR E “ы d mg rep PADO zor 
Å NC КЛ Се Mn Bi he ч AR DINE y ry D 4 FE Ya " DA m Ke gl Eben epe Tl Асе ж ape ее) ea 
, ç А NK IV К 8 ИН" 84 ARE b s LINT ад е бъ est "TT. OH IRI AG a oe. Tor 
E o v. | ° L NG т. i 1 heset UL ЖАҚ N PP КАЛГ аа! TT БИ Hates Y cov 9922464 f Ag Urin енси 
2 n DENN "I ii ы Pee ЕТ] D E EL + 05) ae! LI] M + Мч ыҚ, 724 [> DPF уе TP 4 P p 
А HV A. MT ад Ж КЫРТ, ЖАУ ЖОН 21540 MN e SABE r1 Fio EY TE M DL EUR TD 
ы y ы А Й Ld 1 2 ` 4 d А ЖТТ. in 722% Ld b y Ka А Bef ы 2%% iin rri e am оса 
“© б % “М. wt to Ly Аге. HUU "ж LE LUC DU P АЛ ыы p ым ee Poele RINT "rd polea м 
" : g Н ЖЕК, : ۰ وا‎ > м ууа ) БҮ p 5. d LADE M 01: IE РТА 4 E E аы ke poro — 
) Ü 3° ; ү P =i 42 (257? Мы P bes n ар ost i Jw vnu s rens dte + ДЕ м DE RR 5 A TET э уну IET Ойтал улы Кару балд 
E ее DT JE ПШ ‘ea ы” i F P LET 0 4 Path 4 12 Аг 474 . E rmm Гр б (bv e "nm ee 
' ^ E а ч ЖТ. i IS dg. T = ma; Р 7. Àh ғ)! A a k М4 ты д ые аж, Вор. 22%. 
МЫЛ VAT. Lv di ЛУ Banter rit ти; 21 59 NON TE ETE EFS ERO 
Å 111% UNT Bow y eo ON EU t RO n TTE er RE: 62% IDE: и POE Een RE ten 
" TI E dcum Nap tae te E PERKEN TLEAN RLA SRR S ЫА A eiT T 
Ln - t LU We P Ade a PT А A j ЕД “ea n E io dU HET ft) be 5-45 ALI (7-е оса "Эв be ts СТЕ nz He LEM u æ 6274224 de oie Ws cae ак энин ы — 
d * r 5 7 А + £ ^ T ж 2A 2 М f; ^4 necs M ain m "Um A 511 M жиы iU C Te E E EYN UE 1 E ees 
À uL 1 М "ат V ári, x "Mere 7 m $ y^. y $2: A Уз ) Um PTI EE u ICD dep vie 
Е Е 3 , қ i $ NE шк. R: РЕ, p MP E КОДЕК ai ҺЕ; / d TTD AP OE Jy ei. da йүз; Farið qe hen a NET АІ. o VAT 
3 d ' Á Г кыш Ë ð 3X J t Weis «id ate PITE % E Бб KIT E p^ + Sø pm: ; 
J ж. ei ^ T Lr "YO TR P: В 4. T HS ЕС us mH nr A" "rm MOTA ae sA %,» v Y werde mt C aec 
£ ES. ` Ке M K к eka <i t "n mu TTA ЖАТЫ i} Ж ði ИНН "2 mom U (5 аа Ну dd dt СҢ > ET 
E d dut ы; И Еч 1.3 š Fg V AM vh ite и Mo nr D RM ja ñ ие: Ar LU ief jury rg ra фа pauore eme 
: AA 7 Mak MART TU 2 о s EE a = Lm pm 
о. ‚ ы 0 ч > РАТ ҒЫ TM A ut ; 
i 2 Е ; 4 fe f P Ы IX X. pe E ha M K V Ab th lee ry Hi HN 4 мүге. "epe ë N ON 4 Tp "i 4 Y ag ri я " e Ma ie 656% ERES 
0 1.44 T 7 QUE б NE uw JO UR “мы М. PORE “ыз 5224544 X o هاا‎ "E жү M MI TE j^ IUD те eke " "e od Ағ Tete гс ж 
по E E , ir > s ее o Dee m a a Á d N ac A dhe er "OE eB p s£ t im NO "> SA ME: Ane ú 
E "a n V С А1: ' ^i . 3 crs АДА, JE 4, SL [LE Уой EA мм: 1. PAA Pe ау Nor S Аы КЗ K 
А А = å T : “pad „о % 2 Н $ DOM ۳ Y Я E^ ААҚ e рн y’ L H mh A A Y MA E <a MM ГЕ sia KAJ ree 
: АЕ ИТД i АС кг УС ДУ Dope SERERE КАШИ ud. ATS МЕ rd EIS Ee. 
T Г ds “ТЫҚ i MU "yu Id Л m А. БИП PN) B E Т а Карч УЕ CO 
x de a ТАП “Ж. D K A n к ЖЫ " MILLS HL D PM ur разавай, А зм O зед E y Иа 
AU (SPEM t BAR Dd AY Á An T ERE SE ROUTE E СЕЕ 
۰'9 " 3 (1 9 " Ё UA Жу ETO y 187% б hme Pa ¥ A ET TM H t И UD aug! ” ШЫ Н D F ~ E АТТА n А 
А py ШИЛ fs D S V 4 f ми Á a n " ож $e Алы Ў A a 0 аб gx: A dvi ea У Әр м ПОР) РЕ И Ы (nt 
Á U ei ‹ e NT vel as e SA CALE ТИНЕ CE T 3S a t fati ANE 125 Li d aig La Mitch in dh Te аја P Marks мы 
Ú ^ N " Р ЖҮН” D ы е zm PT ` LI. MI a Ui K 9 Fr АҚҚ желу . mx æa tg Laba ra ae X 
у 1 ч ‚ j An UP BT Жу "vd LM ИЯ IO AE TEL T 417% CREMA US reer Ы nk Ok кетке zet Е 
К ; л m. X p "m TR Dd Pen 75; EA nw Ж ТЕЛДЕ D x eh жігі 5%.” ч DOP ats "i PALA 41 Deg abst ts ы PX ré X cm— - NY? 
А - ь.. ; 4 Г we ҰЗ m. А: N А Ur d OR dH. E e uy: E ACE m x ge n AL eie .. EERE meos = === 
Ы "i, В 57 TT A H iss fiim TK TR 97 Наа H 0 مو‎ 4-1 Аы. ay چ اتا‎ ind 
" D å ИКУ i 3 2: n D FCD (ipo yay rn o EU A [Ls АЯ: a [E TD <... [ кі 58 и 
о, D 2S КЕН y WEEN ák Agt Ly, ҚАТ кың Rn? HE FA = hat aay LR M 
i VE ut d f mM و‎ IPE Domi 1% mt ig Ho 2 s ОЗ? krat HECER 
e £ 2 e SA. . ҢЫ T UT eios 
E vx ie T KUNNAST HU TIR maf, H nt ED e er 
£ ы М И HTETI + جیا‎ Н КТ! И 
FERE vr SEES, 7 finn AE "e On z 
` ТАН T " Mh A U hær par HEN Lt Prot. Mora 
HH "d FEN "Á ЖР к=: EM, FØRT н Lu аз CM vs E IS M aim E. 
Ë dcr le УДА? Ку КЫП DR АКТУ БАРЕ ERG 
; 41/43. CER, prone ; PAR МШЕ VE Pe aks y "E We 43 Aia rk Ta POP зае 
М IST 219 TET 4H р Ч ж ers БАШ, НЕТИ ч pea CU EE E 
, ЖК”, 0 Tt nt Ni r ами RC ДЕР FARA RS tend en 
Š ` ' U { қ eu Мі ris; iin 1% n 52315) AME рага. 
' = ZI 4 A T 
E Ë 3 че ть 
' U 
^ А ? H elg rH 
X К 
Жы; d ea а AA ах” L1 
) э m IIT i rt Cq n ret 
| ыр em 42%: EC Ba lr hi 
| ЧИЛИП. ре ME Se 
NH ut P i, Ен К, Таса 
ie P Ққ ғ bør 
- Ó ee erie 
E de E 
hy e 
vk 
471 аг Æi En Z: E 
iz i VENUE i.e 2 ЖҰ P 
ы x 1 4 Á RE Елі CEE ry E "S 
xn ÓT зн e mic 
; ар, 2 лен E? denn ME RENS АА, 2 
А . TG LER S Ju a eri . 3 а et jc. Wd 
| : b e ens TEN SA L E dto 
i, И E ! J m PES ur Е? 2 Sa Au 
ў i B tt 
1 Ime 7. mE 
dni ni А 
* 
LÀ . 
4 
А 7 
^ š Ø 
pi 
1 j Jeri 
vr 224 
` FOA 
? ы К HE ud 
4 2 
ы E Р +P MEI fame 
AY es 
CL, t en с ae TE 
`. Так ES 5 = 
, nf RON 1: 6454 
4 E Í ia d 
" t Hát 
н 57 
" Н 
p 
p A 
i эю d i 2 
с E 3 Meroen 
ef es Eon 
ا‎ аны Fp а = бач 
Py gon. TE F E E bof ТЕ 
| Ж Ы rU be Wes Ді TÁ MER ч e 
1 i M M d n E D^ 4.2 M tee hy 
Ð ы 9 so ge Тыр por қ 
; прате ла Ü s. ы 
G р 
ЖШ فو‎ A e a E حر‎ ‘Saipan цр 
5 ee} MURS LAMP Үк E mE Р تدوخ‎ ax „кабышы h + 
. s 1 s " x А ra n 
eA = Ex ts 3 SEE us сеня an کج‎ 
`= ie gu: SEL MR Ы AE as „ыд ig ет 2-2 а рЫ E у-у к n a 
v ate ede eT г" dra 424) j = а жа y en ot d кер? Vols “a 
Ж 24) Жы Vim е! rd ad 724 kad ize р, э Б = 
; iae Te pr Pee сатр 4 DA d Hed ll ç: Е а Са UR AUT pent 
М Ao La E q 5514 2 us V "š ag pfe SIT Гіз L 
yr А7 x М "5 4 м “+ Fay ° i fe 2 AE SU. Ets * andi Í e L. 
1 М е a реаты 
j > ШШ, jo M Splints PAN Sida en ед AA e pub erg ici WII = eae 
| ТЕЗ rec pn М H E š w e a 
| : banks SPON da wt jig TA - P eni VA E TTA А | 
А r 1 K d 13 na? Or Nej Jew ` Ж] h НЕ, a bw t. asks À er c H - E رودو‎ » á оа 
Я ED йт ^ E E. 3 4 UR A KE Satis sit FEE 3 TON paren 
J S. ر‎ r ығ $4 4” 16 19:47 IAEA e He 54 x doe ie ЫР Неа Т» ЖЕ ж. b cues 
a ^ tit ROM | ДЕ KR Ø d d. 5 4 TM sb ЕР M ela ted a rox a om petii caa iit 
Ж ° ТІЛЕГІ TA wf vh ml | 2 #\ А TU i % ма ree p иө sð r E 
= r Р. re an H "A M 11 a Ми Muy be derd СА 
„ЖЕ “ыт, ‚ iet. 
, Én a ! ET An 5453 duy rele N М0 ШИ r қ, š ЖУУ BTE fu Mc B ТУ, дігі 188 ab, н у nod ar o rra oe me arden 2 
* é Р > le „ки Ж. e? : Ил OA h ax "o m AUT Р е ("ns OR M ide dae АЗ Betta an? E AE AEA 
ЖҮ ы "d b 5. k $ X ус M Í pt E: LIRE ДАП "t ы ^ k ұғ d еч PS “ЖУ ЛАЛ э AA LE Ke E pears CX re ا‎ Ste d 
^ ; P t E s / م«‎ p ek Я M 0 ЫШ + Ч 5 4 . ê 
қ : g a Ш В 1 к t k ч у е ұй hu Үй: Ki e A^ 1D 2%; 2 S j4 BA 4 қал) Сы iiri tl nk y PM e ed EH ect ie 
ы > ы ۰ ч I t i r vw U T 007 + 2 
, < s D : АМ; 26 ДУ ы d | Ü ТКО E е “ SE ТМ - Vr ы Я k жатық. ES S <; " íM Hond ger Eie ры p q> 
И 2 A xL HED DN ort Eo xe Д | kvi лы hs Sai TUA Куул e А я hed n cr. 
ME e KEN “ЕРА MO cor M MS. MOREM RUE pe uu S s 
Á EC v 3 Халы E ac р 2" à „Ж Ke, rat D. М að Ý в ТАИ e DAN ДЫ RUNE yv. eU. _ er а, с. t ace 
" a 4 u? е , 97249 iip ses n i "Pa a rants “V ía. d E ж varn toy AY E IR "PA a e EET 
ы ' : ' З А pub. t ra #26 P ef P oy £ AN eg. ' t TY Has D ссе Жәй. 
АР . Ж. ж... СЕ ^ d Ux. álin Lala k d ah 4 34$ 4; RAS ras AEN Dare [au Ka a ^ مه‎ dr ane agas la hese ped ashes lyase y а 
' H { " қ jw & - á UL "e» ZUR ID my ПЕТЕ e ? ИШЕ н e EA 5252 E ap as Apgar Teka ht іі у алура араа en 
Ë 3 E BLUR КЕТЕ ЖА PAN aw Wa RA beat: HSS nets EAS Ee A үе ben 
ы + ГЫ ais ' '. L t, # A ды 19%, Ep. ik AL (А 23 Tm Á Er ЫЫ) y y ard بی‎ bd aen m NY e he EC E t oa 
Я 195 5A e 16 d - $t ae m^ EE ALS Pis Mi pee оа TE р rc Mrd AM ен estar heh pees an = pc a 
ү! г А А "uc wn : A é è p D per d p e q Ayla (ААА) У Ду, AMA R s р ر‎ ao y Rd d e r= renare 
5 sa M D Р ч or ES 4 МА. Т "E "ET H [TM * Pn. MITES Used 2,7% 011 ыы ө ы) `: en 5 A c 54 rø 
j ó 4 M ы DUNT CAU 2,4 P hv yon ae tr AE d ым m 4-0 A 574 e abes XN bb ЭЯ, г Die RA tern FRA pap o 
k DEN. 1 gl us ° ! ЗАЧ ы ау iade 1: bg sq xr 117 #' y» еа у 
4 «y hae <, W ne E fave Me Ұмы; P а 0 “ызы T ds t. ep ary e oh a, } NA e E 
3 ы e P d L r ie. yt B eo x 4 к ы б NL $ UN кт y^ PN ۹ x ALS PEN .” 2 vA da: x^ $ ұз 4: E еі ке құма tyd гы STR Ur A TURA alg tre arme. 
{i У dE ONES ES 2 0 A ^ M Y: RC (5552 з he X €x n A sr MR 45: beady жұты; 2 A A cede ета 
b А р و دم‎ с 1 э П x ? E A = A *%®* *.. SKI мА W^ Г H B 
4 | | ИСА a, н ЛА. Ë: E kra ee eri Ар Aqa S; i ca Pian) ERN En en - 
| i dira ce d La алд MC a Ri A ya MAN TENE m 844 қарда 4и А5. Аел Ад аб R ы 
у К J 0 A y sne ы УУУ “ns а A, Fob nd ылы Pp: FØR аш УО ктр" Ый Аай аА etapa en :' 
Я ' А "e Wer O ПИ kah RI AA d pn RA DR IA A k АЗДАР Т мыр ys us TSS rd pre 
( H ы ы Ж = r Т күчө Ф УА” I ra Y YZ ч Daum e esa 
| L , əс е à pee Ps Sat er KR Аг ОН ч, ч ore У каруы e s. MEN yea ace RES 5» Е 4 ج‎ 
А "^, d 6 vn 7.4 € туча M Же, ч Е " 
Ly d À ү К 1 7747,6 RP EX 300379 5 йг! кала ақан DLE MEO e еа 24-2 ^ 
Á р П 4 Гн 440% Жүз و ام ر‎ 2? IT 44. е 09 6 Jara mi Es ane ех SA o gk oo e 
P € г í V W! 244 4 47 Pre 4 ES ча Н HOT re s a 
| P Ü А ољ А Аг AUS Sap ferie MT p P brt /* oj E 
M z | s با‎ ike Ji^ Vra Жр, " z's v 4 x к TN rey de D t^ Je АН V a ee 
5 Ы k ^ q Ч ^ ге 1% ees ør "at. сенщ Ж 7% eio ord (od Adan x) Mri ا‎ 20 -— > 
1 E T! RES RE 2 
| | | | P De RET EE 
. А HEC „* ж, er я vinh LI 
4 : "D s рау қаң 4 Aet 
M А ٤ 1 Mo p i A. hey Coat а PLE 
{ M . e 454 Pests vern ң 2,54 
P М: in еби» Б “ Ық 
i SE LAS беле YS Nees 
a ә = 3 F Ы 
| | : қы Айы ы Ea 
и a g ( by? & n EXIT 
ED i». 
> L] 
k ^ t 
Р қ 
¢ Æ. 
LI 
Ы . 
H e = A L 
› 
1 
, 
M ۰ 
" 
ca 
2 


1 en 
سويت‎ 
^g 
tcv t 
E T 
4244312 
lova har ey 
اوا‎ 
E ree piel Ei Mei. 
4n мач 
^ DR 
ЫС қ pO 
m t bet д ay n 
að Малы 
rah BAR 
à 4 
UT nyt RPE o ЖА 
: k uk 
J 44, р a 
a I + ' Mas 
k > ыла 1 
`. 1] p x E 
Й Ы a 
22 


S ee =< 

Senet te 

4 Å 

Rose Some 5 аера 
pal 

QI LP 

лче T SSeS 

ае [a 

N БЫУ, Ak 

Ане: сара 

e n m pret 

EN PNR 

he af US ру 

wx Fal otn ande le Hr 

TC 425 




















be agate an 
м - verie 
д Деме sv 
ард! эы Б; p A As P 
Aries 29255353 аА a EA LUE DM 
* Ar Nee p^ hA EE кла Ал. ере А E LEE CIAR کاچ‎ 
< ran dir i e im, Saan: 
Od "HET oA sm ا چ‎ be: EEO + р 
662 жар “4 ЕЕ а xls 25 Shree 59% 4 En s RES 
r0 furl An усе ^ UI ONU t EON t eiue Ie 2x — 
RR NT М NOM S D$ COM А, AP ақы» қалі ае Par A t а) чә ҹа pe 
v^ ate aha. ң e оаа ао ра 1 TO РА ا‎ VA, Ta < «4 5 A ا‎ "imd t 
E E en р hA ke ұйы ا‎ an ep A AR ГҮЛ hh ad Cees 251 Ча үйү ме 
ر ور‎ AEN ыы» M. x Vaska 54) {ус з “АЛ, jeu. yy EM tu з-д) = қы хош А рау 
T dur ODES Jo ME AL БТ 4-4” T: VØ R UE Sor уллу ШУН A ажы Hee Te ete 
s Н у ҚАРА K фм uh ENT. мады E Д vbt raa DÀ xA x A tr M | 
nc uw 1059.54 m ê A OR V DREN "Ф: ғақ: JU S Ф Мз EE рр аф узды 
- "mt i | ی جا و لی‎ Ly: ein ғы, m EA а уа SHAR ` E ERE: Ка аьаа 
> ы HA E 97,4 у 4 AP he VM Rhee Hyer ng i rs x, Exe. AA gask وا و‎ ыла ааа, 
4 КҮ Ë % M ы ГА „к ЫДЫ, член: A One pret ae n rn p 
; 54 Балы 3 Ну ioe vb “aps PAGOS ir one e ec wai am کو‎ Berdien as 
a M APR, ш ^ y а 67 vy à ر‎ < % кад Tí ti vehi t$ *h i m = "a A 
کرت‎ SKD Jr RAS. MT " А; e cond v ъч, Ak E >; RA 5,42% ertt 
қ te гө A ONE QE ЖЕК ГАП SA S rhet Мыл Ree = 
dos ° A Es a dne X Akis T би ERE S eR ua SE me реу = a TS 
i H POR UD чар ЗА, 4; К E MUN PORE ГАЛЫ je зды в хы æ E 
. e «* 4 ЖУГУ L^ ерге YE 4” Т T n9, x Eo Жас че 4 
% d £: 25 ki қатар b LX ыы X: чы DANE vm dp aep oh Bee CA ee reme 
М * P а ; hn E Wh > fe ENS NA 4R UU. Din SEES CE fail eer оар Mi 
š pucr Ds М 2 “ы 
` i А е f К а ДАЛЕ Н Н а 3 ft = V RAR 23^ SS Sich лраф: PO UA 
> IR ATA w М, جیا‎ IND we vi MUN Қоры» Сы ауы; 
> M < UM RE. RE ME v AI e 
jl ? З Ð e UN Р LENER KA IIE E AA 
" / sað Á М М р г 
š x j c u VC A yt 4 "s AR gu > е 4 e aaa EPA Lad E 
Б А {у th А tú A: Ae ei te epe 
۰ Y е. veg Sa SR EU: Bae eee te pens ee. 
i t TP қылын Ты, at, 592 e p Pe 2; VOD SUD 
` E ` : ^ Pp ы 44 3 Жа. 4% Алы А. 3 470 ра re ر‎ i p. 
| | R E DE. KLE kar O в 04 
У I i D à ` A e ыз, "t y» T ИГ Бады? Л 
к, 5 LÀ E ы, K 3 4 ^ PIC 
= 127: ; | ^ $e el ME 4 
M J * J 
4 А OE! 
b] М . t, > M ТР” 
` 1 
; z " $ 
£ ^ б Ы " 
% 
LI [] Ж. 
` L4 Ы 
А А Á 
s А R A 
LI LI 
3 4 L 














NAVAL POSTGRADUATE SCHOOL 
Monterey , California 











fN fasi سے پک‎ 
L ` 


\ 


wa. 


Structural Considerations for Aircraft 
Payload Modification- 
P-3C Zero Fuel Weight Increase 


by 


Steven D. Culpepper 
March 1991 


Thesis Advisor: Edward M. Wu 





Approved for public release; distribution is unlimited. 





Jnclassified 
security Classification of this page 


REPORT DOCUMENTATION PAGE 


a Report Security Classification Unclassified 1b Restrictive Markings 
a Security Classification Authority 3 Distribution Availability of Report 
b Declassification/Downgrading Schedule Approved for public release; distribution is unlimited. 


5 Monitoring Organization Report Number(s) 
7a Name of Monitoring Organization 

Naval Postgraduate School 

7b Address (city, state, and ZIP code) 
Monterey, CA 93943-5000 


9 Procurement Instrument Identification Number 


Performing Organization Report Number(s) 
a Name of Performing Organization 6b Office Symbol 
Naval Postgraduate School (If Applicable) 67 
c Address (city, state, and ZIP code) 
Aonterey, CA 93943-5000 


a Name of Funding/Sponsonng Organization 








8b Office Symbol 
(If Applicable) 















c Address (city, state, and ZIP code) 10 Source of Funding Numbers 


1 Title (Include Security Classification’ STRUCTURAL CONSIDERATIONS FOR AIRCRAFT PAYLOAD 


AODIFICATION- P-3H ZERO FUEL WEIGHT INCREASE 


2 Personal Author(s) Steven D. Culpepper,LT USN 
3a Type of Report 13b Time Covered 14 Date of Report (year, month,day) 15 Page Count 


Master's Thesis From To March 1991 77 
6 Supplementary Notation The views expressed in this thesis are those of the author and do not reflect the official 


olicy or position of the Department of Defense or the U.S. Government. 
7 Cosati Codes 18 Subject Terms (continue on reverse if necessary and identify by block number) 
ield | Group Subgroup Aircraft structures, Center section wing box, Euler-Bernoulli, Finite Element model, Structural analysis 


9 Abstract (continue on reverse if necessary and identify by block number 

The Navy is considering the feasibility of increasing the patrol aircraft P-3C zero fuel weight enabling avionics and payload 
rowth. This analysis examines the consequences to the structural requirements of the center section wing box. Two solutions to the 
tructures field equations are investigated: a simplified hand solution for preliminary feasibility calculations and a more precise 
olution for design analysis. Together, the solutions provide a necessary check for the results. The simplified solution employs the 
uler-Bernoulli assumption which generates a set of integrals expressed in terms of the assumed displacements. These integrals, when 
ombined with simplified geometric shapes and symmetry, ultimately produce a decoupled matrix solution. The precise solution uses 
PC based finite element method which simultaneously solves the field equations for basic elements to be linked together with the 
ppropriate boundary conditions. For the current 135,000 pound gross weight 1g load condition, the internal stresses calculated by 
nite element are in accord with those by simplified hand calculation. Extensions from this modeling will generate design criterion 
or the target 95,000 pound zero fuel weight aircraft as well as alternate flight or taxi conditions. 


0 Distribution/Availability of Abstract 21 Abstract Security Classification 

unclassified/unlimited [| same as report ü DTIC users Unclassified 

2a Name of Responsible Individual 22b Telephone (/nclude Area code) 22c Office Symbol 
ward M. Wu (408) 646-3459 67 

D FORM 1473, 84 MAR 83 APR edition may be used until exhausted security classification of this page 


All other editions are obsolete Unclassified 


Approved for public release; distribution is unlimited. 
Structural Considerations For Aircraft Payload Modification- 
P-3H Zero Fuel Weight Increase. 
by 
Steven Drew Culpepper 


Lieutenant, United States Navy 
B.S., United States Naval Academy, 1983 


Submitted in partial fulfillment of the requirements for the degre of 
MASTER OF SCIENCE IN AERONAUTICAL ENGINEERING 
from the 


NAVAL POSTGRADUATE SCHOOL 
March 1991 


ABSTRACT 


The Navy is considering the feasibility of increasing the patrol aircraft 
P-3C zero fuel weight enabling avionics and payload growth. This analysis 
examines the consequences to the structural requirements of the center 
section wing box. Two solutions to the structures field equations are 
investigated: a simplified hand solution for preliminary feasibility 
calculations and a more precise solution for design analysis. Together, the 
solutions provide a necessary check for the results. The simplified solution 
employs the Euler-Bernoulli assumption which generates a set of integrals 
expressed in terms of the assumed displacements. These integrals, when 
combined with simplified geometric shapes and symmetry, ultimately 
produce a decoupled matrix solution. The precise solution uses a PC based 
finite element method which simultaneously solves the field equations for 
basic elements to be linked together with the appropriate boundary 
conditions. For the current 135,000 pound gross weight 1g load condition, 
the internal stresses calculated by finite element are in accord with those by 
simplified hand calculation. Extensions from this modeling will generate 
design criterion for the target 95,000 pound zero fuel weight aircraft as 
well as alternate flight or taxi conditions. 
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I. INTRODUCTION 


The Navy is considering the feasibility of increasing the patrol aircraft 
P-3C zero fuel weight enabling avionics and payload growth. An 18,000 
pounds payload increase from the current zero fuel weight of 77,200 
pounds will undoubtably require airframe structural modifications. This 
payload increase to the present configuration will accommodate system 
upgrades into the next century. 

This analysis will provide the tools to determine whether the existing 
center section wing box can accommodate the additional loads or if a 
modification is required. Modification options ranging from variations of 
the existing geometry to materials substitution such as composites will be 
discussed. 

A Structural analysis, which provides a basis for strength and stability 
assessment, begins by modeling the geometry and boundary conditions of 
the center section wing box for a given flight condition. An idealized 
model constructed from rods and sheets replaces the center section wing 
box. This model approximates the geometry of the structure while 
providing a necessary simplification that facilitates the mathematical 
solution to the solid structures field equations. The wing boundary 
conditions are provided by a contractor while a component inventory 
method (Appendix A) is used to construct the fuselage loads. 

Two solutions to the field equations provide a necessary check on the 
results. The Euler-Bernoulli assumption uncouples the field equations and 
leads to an integral representation of stresses in terms of the boundary 
condition resultant loads applied at the centroid. The idealized model 
geometry (Figure 1) allows the boundary value problem to be integrable 
thereby yielding a closed form of the solution. From this closed form, a 
solution to uniformly distributed (centroid) forces and moments is 
practical. The second solution uses the Finite Element Method (FEM) to 
simultaneously satisfy the field equations for basic elements or geometries. 
These elements are combined to form the structure of interest. This 
numerically based solution method allows a solution of more realistically 
distributed forces and moments as well as refined geometric 
configurations. The problem formulation and computational 
implementations for the current application are verified for the idealized 


model using closed-form hand calculated results. The Saint Venant 
principle enables the resultant boundary condition loads to be modeled at 
the centroid thereby verifying the (lower bound) internal stresses predicted 
by the hand calculations for the static case. Upon verification of the hand 
calculated results, the finite element model may be used to examine various 
realistic inertial loadings within the center section wing box for any given 
flight condition in the operational envelope. 


FS 695 
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Figure 1. Idealization of the Center Section Wing Box: rod and 
sheet structure replace actual P-3 geometry to facilitate closed form of 
mathematical solution. 


Validation of the boundary condition stresses generates confidence in 
the solution of the center section wing box for various inertial loadings. 


New boundary loads for the proposed flight conditions of interest will help 
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identify critical stress and deformation locations within the structure. The 
finite element software lends itself to parametric studies of alternative 
structural configurations in addition to exploring the effects of increasing 
the load requirements on the airframe from, for example, 3 to 3.5 g's. 
Other materials including composites may be examined as replacement 
options. New configurations can be readily analyzed by fine tuning the 
boundary load conditions. These benefits lead one to the utilities of this 
structural analysis. It provides a readily accessible feasibility check of the 
existing structure and a means to design alternatives for the center section 
wing box. The analytically generated data may form the basis for Request 
For Proposals (RFP) which may include innovative designs (ie. composites) 
and structural modifications. 


II. Problem Definition 


The future P-3 version H will accommodate an increase in the zero fuel 
weight from 77,200 to 95,000 LBS. The weight added will occur within 
the fuselage section of the aircraft. The proposed solutions to this 
structural problem include replacing the entire wing, and replacing all or 
strengthening weak members of the center section wing box. Each of these 
solutions requires an in-depth knowledge of the structural limitations of the 
current center section wing box. Specifically, the question remains 
whether the modified structure will support the added stress and strain 
given the increase in zero fuel weight. 


A. BOUNDARY CONDITIONS 

The worst case conditions for the center section at maximum gross 
weight include a 2g taxi bump on the ground and a sustained 3g turn in the 
air. The lift generated by the wings accelerate the 155,000 pound airframe 
with a load factor of 3.0. This lift transmits a moment and a shear into the 
center section at wing station 65 (refer to Figure 1) The forward and aft 
fuselage sections impart moments and shears into the center section 
fuselage stations 571 and 695 respectively. The taxi load is generated by 
the main landing gear impulses subjecting the inboard nacelles to vertical 
accelerations. While this condition delivers a shear, the moment is of 
opposite sign to the aerodynamic load yet it still acts on the center section 
at wing station 65. For purposes of this study, the 1g in-flight loading 
condition was analyzed in detail. Other loading conditions were provided 
to cover the full spectrum of situations. 


B. FUSELAGE LOADS 

The fuselage loads for a 135,000 pound aircraft were constructed using 
the data provided by the current P-3 version C Specification printed in 
1982. The weights of various components (ie. wing, propulsion, body, tail, 
electronics etc.) were broken down into fuselage weight and wing weight 
both inboard and outboard of wing station 65. The fuselage weight was 


then distributed along the airframe based upon the location of the 
components. The shear and moments were built up from this weight 
distribution and tabulated in Appendix A. The loadings along the center 
line of the fuselage were computed from: 


m 
Men) ae (ОЖ) 
1-1 


тп 
Уау w; (2.2) 
i=] 


w; = weight of the 1,,component 
l; = moment arm of the 1,, component 


n = load factor 
m = number of components 


The moment arms were referenced to F.S. 571 for the forward section and 
to F.S. 695 for the aft portion of the fuselage. 


C. WING LOADS 

The wing loads were provided by Aerostructures, Inc. located in 
Arlington, Virginia (Appendix A). Their data was extrapolated from 
values given for the eight loading conditions in Lockheed's Structural Life 
Extension Program report for the P-3C. The bending moment and shear 
were provided for 1g and delta 1g increments along the entire wing 
starting at the center line and continuing out to wing station 584. The load 
accuracy was quoted as between five and ten percent. The moment 
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contribution about the wing My is calculated as the shear multiplied by the 


distance between the lift line and mid-chord (approximately 0.2 chord). 

A first order method to calculate the wing loads requires the lift to act 
through the center of pressure on the mean aerodynamic chord (Figure 2). 
This moment and shear felt at wing station 65 1s countered by the weight of 
the wing and fuel (wet wing). The total moment and shear at W.S. 65 1s 
determined by dividing the wing up into sections and summing the 
incremental shears and moments produced at W.S. 65 by the outboard lift 
distribution. 





Y | WS.65 WS. 254 WS. 571 
mac 


Z 


Figure 2. Two Dimensional Lift Distribution : simplified 
resultants (as seen by center section) shown at mean aerodynamic chord 
(mac). 


D. GROUND TAXI LOADS 

The ground taxi loads were provided by Mr Nam Phan 
(NAVAIRSYSCOM). The maximum gross weight condition occurs when 
the plane rolls over a 2g bump and the wing responds in a flexible manner. 
The gear through the inboard nacelle accelerate vertically upward while 
the fuselage and fuel ladened wing outboard of the nacelle resist the 
motion. The bending moment located at wing station 65 during this 
condition creates tension on the wing top and compression on the bottom. 


E. CONTROL SURFACE LOADS 

A rudder kick produces a torque about the fuselage, a moment about 
the vertical axis and an insignificant y-direction shear. The vertical 
distance from the rudder center of pressure to the longitudinal (x) axis 
multiplied by the lift produced by the rudder due to a deflection into the 
slipstream adequately describes the torque experienced by the fuselage. 
The moment generated by the rudder is a product of the rudder lift times 
the longitudinal distance from F.S. 695 to the rudder center of pressure. 

Elevator deflection incrementally alters the shear and moment 
produced at F.S. 695. The elevator force directly adds to the remainder of 
the shear at F.S. 695. The control force was multiplied by the longitudinal 
distance measured from the elevator aerodynamic center to F.S. 695 to 
produce this moment. The moment was added with all the rest at F.S. 695. 

Ailerons generate incremental amounts of lift in comparison to that 
produced by the remainder of the wing. As such, they can be readily 
incorporated into the calculations. The shear adds directly to that already 
calculated at W.S. 65. The moment results from the lifting force on the 
aileron acting about the moment arm established by the aileron center of 
pressure to W.S. 65. 

This summary of load conditions exhausts the list that need be 
examined for purposes of structural integrity. These conditions should in 
fact be pared down further to a minimum set that can be quickly 
incorporated into the model and determine whether or not the structure has 
been overloaded. The answer to the overload question will determine the 
need for wing box redesign thereby closing the design loop. 


III. GENERAL SOLUTION TO P-3C 


The general engineering solution provides a linear analysis to the 
question raised regarding the P-3H structural response to an increase in 
zero fuel weight. This solution yields an essential comparison for the 
output of the finite element program as an accuracy check. This 
formulation is based on an assumed general form for the displacements 
with the parameters to be determined by the specific boundary conditions. 
This formulation reduced the solution of the field equations to algebraic 
forms (after the necessary integrations) thus bypassing the necessity, of 
solving 18 partial differential equations simultaneously. 


A. SOLUTION USING EULER-BERNOULLI ASSUMPTION 
The assumed displacement method was originally used by Euler and 
Bernoulli over a century ago. Essentially, plane sections are assumed to 
remain plane during application of bending loads. The limitation to the 
theory lies in its application to only small deformations. Strain, defined in 
terms of the partial derivatives of displacement, can be calculated for the 
Structure in terms of the assumed displacements. The constitutive 
relationship for a specific material (i.e. aluminum) converts the strain to 
internal stress. The internai equilibrium equations establish a set of 
equations in terms of integrals for solution in terms of internal stresses. 
These internal stresses are then equilibrated to the boundary tractions. А 
system of integral stress equations related to the known applied boundary 
conditions result. The solution of the integrals in this system of equations 
is further simplified through the idealization of the actual structural 
configuration by interconnecting bars ( to carry normal stress) and sheets 
(to carry shear). The structure in Figure 3 models the center section wing 
box of the P-3C. The internal stresses within the individual members are 
constant permitting their extraction from the integrand thereby facilitating 
the integration, which results in a system of algebraic equations. The 
general solution, now represented in matrix form contains many coupling 
elements. Decoupling is accomplished through the application of symmetry 


to the geometry and using the specific loads that apply for a particular 
flight condition. We note that this simplification is possible only for the 
simplified centroidal (uniformly loaded) boundary conditions. 
Nevertheless, this solution provides a useful estimate of the best case 
expectations (lower bound stress) and it also provides a bench mark 
verification for the finite element method solution which will be used for 
examination of the realistic cases in terms of payload distribution, 
structural modification and structural material substitution. 





X 


Z 


Figure 3 Center Section Wing Box Structure : Each circle 
represents a node of the structure. Shear sheets exist on the vertical x-z 
and y-z planes of each cell. Rods are represented by the lines which 
connect each node. 
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I: 


Summary of Field Equations 
The five governing field equations are as follows: 


Stress Boundary Conditions 
O yx O xy O xz Vx 


T, 
Ту |ғ| бухбуубу | | Уу (3.1) 
T, O rx O zy Oz У, 


Differential Equations of Equilibr m 


00 xx " dO yy z ос, 





+Х =0 
Ox dy dz 
do oo do 
EC 
Ox dy dz 
905 , zy , z 7-0 
Ox dy dz 


C: astitutive Relation 


І-УУ v 0 0 0 
v 1-v v 0 0 0 
| Е v v i-v 0 0 0 
KAV lV) 0 0 0 05У0 0 
| 0 0 0 О 0.5-v 0 
00 0 0 0 0.5-у 
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(3.3) 


Strain Displacement 


(3.4) 


Strain Compatibility 


Den den дЕ 
2 БА + уу шы E 
д | Ey , Erz Exy | _ 20 Ez * 3,2 Әхду 
dz\ ox дхду 


ду dz 
| 2 den OE dey, 























дх dy 92 дуд2 dz dy 
2 2 2 
д 9Еу, _ ОЕ ‚ ОЁлу _ 20 Eyy д Еп, 0 Ex _ 0 Ex 
dy\ dx dy dz дхд2 s ДА дхд2 


Тһе stress boundary conditions, Equation 3.1, relate the external 
traction (B.C.) to the internal stresses in the first internal layer. The 
second system of equations (Equation. 3.2) relates the first internal layer to 
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the adjacent layer by balance of internal forces; Hence it is known as the 
differential equation of equilibrium. Equations 3.3, the constitutive 
relations, convert strain uniquely into stress. Each material exhibits its 
own particular constitutive relationship. Aluminum behaves in an 
isotropic, Hookean manner. The strain displacement equations (Equations 
3.4), relate the strain to the partial differential of the displacement. 
Finally, Equation 3.5 ensures that the second derivative of the strain is 
continuous or that the material will not separate anywhere within the 
interior. Any solution to a solid structural problem must satisfy these five 
sets of equations. 

One approach to the problem requires the simultaneov soi on of the 
eighteen coupled, partial differential equations (PDE) described in 
Equations 3.1 through 3.5. Each specific geometry of the structure and 
boundary condition requires a different specific solution. This method 
requires considerable mathematical complexity, and in fact, is frequently 
intractable for realistic structural configurations. The Finite Element 
Method (FEM) provides a general solution to a sub-geometry or element. 
The computer uses its high speed to solve the general case for each 
element. These elements may then be assembled to approximate any 
specific configuration of interest. 

An approximate approach to this structural problem is based upon the 
Euler-Bernoulli assumption. This formulation assumes the functional form 
of the displacement, 1.e., that plane sections remain plane in the structure 
while deforming under normal and bending loads. The assumption 
guarantees that the compatibility relations, Equations 3.5, are satisfied. 
Deformations of the structure must be small to ensure that the small angle 
approximation for the tangent applies (ie. less than 15 degrees of shear 
deformation yields an error of less than 3%). When the form of the 
displacement function is assumed, the solution of the eighteen PDE's 
reduces to a set of algebraic equations of differentials. In this case the 
displacement is of an assumed functional form with parameters to be 
determined by each specific boundary condition (B.C.) applied to the 
model. 

This approximate approach begins by differentiating the assumed 
displacements using Equation 3.4. 
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у(х,у,2) ту(-0,(у)2-0Хух (3.6) 





а а 
- 9v . Wo d X Ө; (3.7) 


The unknowns in the displacement relation become parameters to be 
determined. The strains are then expressed as internal stress using the 
constitutive relationship shown in Equation 3.3. Simplification of the axial 
stress expression in the y direction yields Equation 3.8. 


d 
Oyy =Ee,, = E (4У0 -Z е; 
ду ау 





do 
+x— (3.8) 
dy 


Substituting these internal stress expressions into Equation 3.2 yields a set 
of differential equations. Equilibrating the internal stress to the boundary 
tractions, Equation 3.1, for a specific geometry results in a set of integral 
equations. The integral contains the unknown stress expressed in terms of 
the parameters of the assumed displacement. Equations 3.9-3.11 equate the 
general B.C.'s applied at wing station 65 that have an axial stress 
component. 








d dg 

P„=- Í E( 2-7 + x JdA (3.9) 
dy dy y 
d dg dg 

M = | EC 0 , 7+ x 7dA  Q.10) 
dy dy dy 
d do do 

M,s-] EC—-z—'4 x—5xdA (3.11) 





dy dy dy 


A 
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From Equations 3.9-3.11, the model geometry was simplified by replacing 
all of the structural members by sheets that carry shear and rods that 
support axial loads. Through this modeling, the internal parameters 
become constant and can be factored outside of the integral expression. 

















d dg dg 
P,z-Ej-9| Ё. ал+; -*} Paaa-Ej—:| Ехал (312) 
dy E; dy E; dy E; 
A A A 
d dg 2 dg 
M, = و‎ Mo En dA +E, = En dA-E,; а =e dA (3. ! 3) 
dy E; dy E; dy E; 
A A A 
d dg dg 2 
М, = 5 ET dA +E,— E dA E Dx dA (3.14) 
dy E; dy E; dy E; 
A A A 


The resulting equation was integrable since the integrand is solely related 
to the geometry of the cross section. The following modulus of elasticity 


weighted properties are defined: 


A = |P EE 
E; 





x = l A (3.16) 
A li 





pe ТЕЙ ОТК 
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* 2 

basf Ez dA (3.18) 
aum 

І,, = 3 aA (3.19) 
i. 

» 2 

l,2-.| Ex aa (320) 
LI 


Using these definitions, Equations 3.12-3.14 transform into: 


dvo, * dg — dg t 
Ру--Ер--А +E, Z - Е г— х (321) 

















dy dy dy 
dvi" do < do < 

M S F sn p ч и ү (6120) 
dy dy dy 
dv = do < do < 

E cup E TT OSA 
dy dy dy 


Solving for the derivatives of the displacement functions and inserting into 
Equation 3.8: 





P МІ, +М,І M.L, +M 
ES ‚ E‚MdmtMde, E‚MwztMde, (324) 
E E; * * „2 E; * ж .2 
ОШ же Ш>) Gd de) 


The differential equations have thus been transformed into a set of coupled 
algebraic equations. The shear stresses are solved for in a similar manner 
and attached in Appendix B. The general shear term 1s: 
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0 (-m,- V )I,, * (-m, * VI Q 
Oxy = Oy, - + صا‎ рн (3.25) 
(Iy lz- Iy) 
S1 
i E 
О, = —y tds (3.26) 
1 
0 
S1 
* E | 
Qe — ‘ds (3.27) 
E, 
0 
The resulting stresses for the general case are: 
Б. MAy +M M Iz +M 
om _Ё _ E May tM yy Ep Mra tMdy (3.28) 
Ei E ж ж ж 2 E; ж ж ж 2 
(IyyIz - 15) (Iy lz- ly) 
P MJ,, -M,I M.I, +M 
o, BT. NM de, Ep MdatMdn, 399) 
yy É 2 2 
E1 E, * * * E; * * * 
Med Cl dE ede) 
P M.I „+M MI, +MI 
o, EV: E (330) 
Ed E; * * * 2 E; * * * 2 
C Ij Ix - Ley ) (Ia) 
-m, € V JI, + (-m, - V )I 
oncom _q(0) _ RS (3.31) 
t = 


(hy I 
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-m,- VL, , + (-m,+ V,)I 
R — азу 
t е. t 
ЕЕ 12) 
E _q(0) _ (-m, +V Jl, + C-m;- Vgl, О, 
Oy = Oy EE en oe (335) 


A 


This system of equations can be arranged in matrix form: 


о] = [к] [F] (3.34) 


The nine generalized forces consisting of three forces, shears and moments 
are weighted by the k matrix to yield nine stresses. The terms that make 
up the non-zero elements of the matrix are listed below. The generalized k 


matrix follows these terms. 





E, - I. oie 
ky) =k52=ko3= - » Kig= у = 
E A E; * * „2 
! (Iy1z - 17) 
E. -ylI, + -I Q 
Eee e] ; pond ec оъ ету 
E, * * „2 * * „2 t 
(Iy,Iz - L7 ) (ly Iz - 14) 
Г Q. г О, 
k= kg = , = 
LE) (IyyI7- Ig ) 


1j 


= = = E 
K 36 K% * * 2 ( <97 E. * ж ж 2 | 
(Ij Iz-1I,) Cl at) 
M E х1, -21,, x = О, 
Kss ;l» ke-kg- —: 
Ch 1208) GLa x2) 
Di * I. I І, 
Коб = Кв6 = = 5 = ; к-г ee 1 
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Further simplification is possible by applying the known loading 
conditions for a given boundary. For the y face at wing station 65, the sole 
axial load is a moment about the x-axis produced by the lift acting on the 
wing. Accordingly, Equation 3.24 reduces to: 


M,I М, 
O yy БЕШЕНЕ Ее У ELO UNE MEA (3.35) 
E, * * x: 2 E; ж ж Ф 2 
(Ix Iz- Íx (Ixx Iz- L. ) 


The profile geometry of the wing box is not quite symmetric but will 
be assumed so for the purpose of this calculation. Therefore the cross 
product of inertia terms vanish and Equation 3.35 reduces to: 





(3.36) 
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Similar arguments hold for the other stresses in Equations 3.28-3.33. 
These last two simplifications lead to an uncoupling of the algebraic 
expressions in the matrix as shown below: 
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Each stress relates to an applied boundary condition weighted by a k 
factor determined from the geometry of the cross section. The in-flight 1g 
condition generated the first load condition for the model. My represents 


the moment generated by the fuselage forward of F.S. 571 and aft of F.S. 
695. For reasons of stability and equilibrium, the fore and aft moments 
аге equal. The My moment is produced by the lift acting on the wing 


through its center of pressure, the quarter chord position on the mean 
aerodynamic chord. Two sets of vertical shear exist. Established by the 
lift on the wing, shear passes into the center section at wing station 65 
acting in the negative z direction toward the top of the * selage. The 
second shear results from the fuselage forward and aft of the nter section 
wing box. This shear acts in the positive z direction at F.S. 571 and F.S. 
695. P, accounts for the mass contained in the center section itself. 


It should be noted that this solution, even with the Euler-Bernoulli 
assumption is tractable only for the hypothetically uniform loading where 
the resultant forces and moments act on the centroid of the cross-section. 
In fact, there exist no structural component at the centroid to carry the 
applied loads. Nevertheless, because of the mathematical tractability of this 
hypothetical case, it provides a bench mark for verifying the subsequent 
finite element method procedures. This idealized centroidal loading also 
provides the estimation of the best case analysis, i.e., the lowest possible 
internal structural stress for a given payload. 
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IV. HAND CALCULATION 


In this chapter, the specific geometric properties of the P-3C aircraft 
are reduced to appropriate forms to be consistent with the Euler-Bernoulli 
formulation derived in Chapter III. 


A. WING MODELING 
Starting with the view from the left wingtip looking toward the 
fuselage, the cross sectional properties are modeled from the Structure 
Wing Fuselage Intersection diagram found in Lockheed Report No. 13102. 
At wing station 65, a double shear splice connects the upper and lower 
wing surfaces to the center section wing box. Tension type splices join the 
beam caps together [ref 4]. At this wing station, the cross section may be 
modeled as a rectangle of thickness t. The outer length runs from Fuselage 
Station F.S. 571 to F.S. 695 or 124 inches. The inner length runs an inch 
shorter. The outer height and inner height are 18 and 17.5 inches 
respectively. The moment of inertia for this geometry is: 


3 3 4 
1 = (bh, -bj) — 5330in (4.1) 


The area is: 
A-4A;-bjh,- bj, 279.5in (4.2) 


The distance from the centroidal axis is: 


h3 = aj I» — 8.19 in (4.3) 
А 


The cross section transformed into rods and sheets (Figure 4). 


(A: 
| 2h3 2 16.38 in 
M 


Figure 4. Wing Modeling: wing station 65 
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Moving towards centerline from W.S. 65, the fuselage must be 
accounted for in the geometry. The curved I-beams that form the main 
frames constitute the overwhelming majority of the additional resistance to 
the bending moment. This addition can be modeled as two rods located 
above the wing box cross section already completed and two rods at the 
same location as the lower part of the wing box (Figure 5).The area (Aj) 


represents one half of the cross sectional area of the seven frames existing 
between F.S. 571 and F.S. 695, А; equals 3.5 in2. The vertical distance 


(h4) from the top of the wing box to the I-beam varies in the real 
structure, but for simplicity will be modeled as 64.4 inches (see fvselage 
modeling). The distance from the wing box to the lower seven fr леѕ 15 
24-2.0 inches. The new moment of inertia term is expressed in Equation 


4.4. 


h4 = 64.4 in 


4 
_ = 2h3 = 16.38 in 
l С { 2.0 in 





Figure 5. Wing Modeling: wing station 54 


= R (A boxZbox + A curveZ curve t 2A g;) 


Z = -333in (4.5) 
SPEI (Atoxt Acuve+ 2A) 


E 2 T 
lxi Ixxpeam * 7. beamÅ beam + Lxxpox +Z box A box + Íxgjaw t Z úlow Á low 


1. 2 04323004 5330 +1730 +0 + 1160 2 40500 in^. (4.5) 
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From Chapter III, the reduced form of the axial stress along the wing сап 


be expressed as: 
|. MyZ 


б 





yy (4.6) 


XX 


The bending moment increases by the shear loads at the W.S. 65 and those 
additional loads along the structure: 


M, (y) = M, (65) eV zy) dy (4.7) 


The last term in Equation 4.7 which requires an explanation is the z 
location determined by the position of the member in question with respect 
to the modulus weighted centroid of the cross section. 

The shear stress follows from the reduced form of Equation 3.26: 


= -V, О, 


xx “ч 


Су; (4.8) 


The moment of inertia term was resolved for the axial stresses апа will be 
reused here. The centroidal term from Equation 3.19: 


CE 
© EA (4.9) 

ГЕ 
As with Equation 4.7, the Zz; term represents the distance between the 
centroid and the ith element in question. Aj represents the area of the ith 
element. The shear load V,(y), increases from the wingtip to W.S. 65 as 
the wing generated lift dominates the weight of the structure, ordnance and 
fuel. Between W.S. 65 and the centerline, the shear is modified by the 
weight of the structure and fuel of the center section. The shears at wing 
stations 65 and O were provided by a private contractor. The data will be 
linearly interpolated between these two points to arrive at a shear value for 
a given y location. 
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Equation 4.8 assumes a centroidal loading which in the case of a wing 
with lift acting through the 30 percent chord line must be modified by the 
superposition of the stress due to the torsional loading. For single cell 
closed sections such as wing station 65, the stress expression due to torque 
1s expressed in Equation 4.9. The enclosed area, torque and sheet thickness 
are included. 


= 


бус 


У (4.10) 


> 


For the case of multi-cell structures refer to Allen and Haisler s text, 
reference 2. The authors provide a method whereby cell deformations in 
an n-cell structure are equated to generate n-l equations. The final 
equation comes from the equilibrium expression in which the moments 
about the longitudinal axis are summed and equated to zero. 


B. FUSELAGE MODELING: 
The axial stresses in the x direction contend with the geometry shown 
in Figure 6. 





Figure 6. Simplified Fuselage Cross Section 


Dissecting Figure 6 further into simpler geometries enables a 
component build-up approach to the section properties. Viewed from the 
nose of the aircraft looking back at F.S. 571, the upper shell and stringers 
of the fuselage may be modeled as a thin sheet (Figure 7). The center piece 
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may be modeled as a rectangle (Figure 8) and the lower I-beams and 
fuselage as two rods. 







Ra =67.75 In. 


O= —л/6 









©= 71/6 < 


TOLIN SRELE MODEL 


Figure 7. Upper Shell and Stringers of the Fuselage 
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= 1 B R, Ө = 7П/6 | K | 
z= J zdA= fr, E (r sin®)r агае = -28.05 in (4.12) 


6»71V6 


I = | саа = Re (z in9) r dr de 4.13 
yYcurve J, 2-7. FØR z-rsnO)rdr (4.13) 


Q=- 


L youve = 1.298 E+05 - 1.120 E+05 + 0.560 E+05 = 7.38 E+04 in’ (4.14) 
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130 їп. 


Figure 8. Center Piece of the Fuselage 


3 3 
[yox = +> (b hi -bjh;)216.0E403 (4.15) 


Ax,-bjh-bd.- (18. (130in)-(16.5inX(126in) 2261 in =8А; (4.16) 


The moments of inertia are computed for each piece about its respective 
center of gravity. At this juncture, the structure converts to rods and 
sheets using the computed areas and moments of inertia to calculate the 
vertical separation (Figure 9). An Ay of five square inches represents the 


area of the I-beams under the wing box and the fuselage shell lumped 
together two inches below the center section wing box. 





Figure 9. Idealized Fuselage Cross Section 
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Zourve= AJ p =32.19 іп (4.17) 
I 
бол = л | = 7.83 in (4.18) 
box 


Solving for the modulus weighted center of gravity for the entire cross 
section enables the section moment of inertia to be calculated. Figure 10 
shows the dimensions of interest for the complete cross section. 


= (A boxZbox + À curveZ curve * 2A zy) 
(Abox+ Acuvet 2Ax) 


=0.65in (4.19) 


T 


64.38 in 


Z structure — 


15.66 in 





Figure 10. Idealized Fuselage Cross Section : with dimensions 


2 2 2 
Lyytota] = Í yycurve + Z` curveÅ curve t Iyypox + Z box Á box + [уу + 2 low Aıow (4.20) 
Lyyrotaj = 7.38 E+04 + 7.68 E+04 + 1.60 E+04 + 1.34 E+04 + 0 + 0.29Е+04 


=1.83 Е+05іп (4.21) 


The simplified x direction axial stress equation follows: 


o. Tt (4.22) 


yy 
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The moment of inertia term has just been solved for the homogeneous case 
in Equation 4.21. The location of the bar in question, referenced to the 
modulus weighted centroid, determines the value of z. The boundary 
conditions at F.S. 571 and 695 were developed in the boundary condition 
section and calculated in Appendix A. To determine the moment at a 
particular fuselage station between F.S. 571 and 695 requires a similar 
calculation based upon this component build up method. 

The shear stress for the face at F.S. 571 simplifies from Equation 3.24: 


(4.23) 


бу, = - 


27,0, 
t 


i yy i 


The terms in Equation 4.23 follow directly from the explanation for 
Equation 4.10 where: 


LE 
Оу = 2, iA (4.24) 
i 1 


The moment of inertia term was previously solved for in Equation 4.20. 
The thickness depends upon the location of interest. The shear force V; 


was calculated for F.S. 571 in Appendix A. All of the 1g boundary 
condition loads are shown in Figure 11. 
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Z 
Ning = 5.97E+06 in-lbs Vwing = 24,400 Ibs 
rear 
Mwd = 3.31 Е+06 in-Ibs Midd F 12,600 Ibs 
Inertia- 2" Vwing LEONE MEI = 27500 Ibs 


Figure 11 Boundary Condition Loads for Center Section Wing 
Box : at 1g. 


In this chapter, the actual P-3C aircraft geometry is reduced to 
idealized rods and rectangular sheets which are statically equivalent to the 
actual load and moment carrying capacity of the aircraft. The actual load 
and moments are reduced to centroidal boundary conditions as shown in 
Figure 11 for the 1 g condition. Stresses calculated from the 1 g condition 
can be simply scaled to any multiple for feasibility studies. 
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V. FINITE ELEMENT METHOD 


The finite element method provides a direct solution to the field 
equations (as described in Chapter IT) without the simplification of Euler- 
Bernoulli hypothesis. For the center section wing box it can provide 
improved accuracy over the hand calculated solution for centroidal loading 
as well as for other more realistic cases which are beyond hand 
calculations. Small or finite elements are solved in the general sense to 
simultaneously satisfy the eighteen field equations. These building blocks 
are subsequently combined and matched to form the structure or geometry 
of interest. Equilibrium must be satisfied at each node where the elements 
connect. The external displacement and force boundary conditions are 
applied to the structure. The solution to this combination of elements and 
boundary conditions satisfies the minimum energy principle. This 
principle guarantees a unique solution which minimizes the strain energy 
preserved in the structure as it deforms under loading. Model One was 
constructed to compare this method of calculation to that accomplished by 
hand calculations with centroidal boundary conditions. (See Appendix D 
for FEM examples) 

Upon verification, subsequent load conditions can be constructed to 
model a flight condition of interest (i.e. 3g, 3.5 g ... etc.). The resulting 
internal stresses will identify iocations which require redesign. Modified 
configurations (1.e. fuel/payload) cause a redistribution of the inertial loads 
within the center section. This model provides the tool from which to 
explore these modifications to determine their effect on the internal stresses 
of the wing box structure. If critical stresses are identified, the model 
lends itself to study the replacement options for the structure. Redesign 
options include additional area/mass added to a component or, on the other 
hand, part replacement with an alternative material (i.e. different 
composites). Both of these options require assessable software 
modifications to the model definition. 

The model to be constructed for finite element analysis is not 
geometrically identical to the actual structure of the P-3C aircraft. A 
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geometrically exact model would be too complex to implement and too 
time consuming to execute; it is appropriate only for the actual design 
stage. À geometrically simplified model retaining the pertinent structural 
features is constructed herein; its application is intended for problem 
identification and assessment of different options of payload, distribution 
structural enhancement and modifications. 

The finite element software used in this investigation is the PAL2 V 4.0 
by MacNeil Swindler Corp. This software can be run on PC class of 
computers, it has provisions for analysis of composite materials and can be 
executed under NASTRAN on mini and mainframe class of computers. 

This chapter describes the geometric definition of the model and the 
software related procedures to assure the appropriateness of the finite 
element method implementations. 


A. PROCEDURES 
1. Geometric Definition 

The nodal points that describe the structure of interest are defined 
Using a rectangular coordinate system (polar and spherical are available). 
Once the nodes are in place and the material specified (i.e. 2014 T-4 
aluminum), the elements that make up the structure are connected to the 
nodal lattice. Rods join two nodes together while quadrilateral plates 
connect four nodes. The circular rods are defined by inner and outer 
diameter. The shear area 1s nothing more than a correction term for the 
parabolic shear distribution across the beam face. The square bar shear 
area equals the cross-sectional area divided by 1.2. Their lengths are 
prescribed by the lattice. The shear plate thickness is specified, while the 
node locations determine the planer dimensions. The relation of this model 
to the actual structure is shown in Figure 1 and Figure 12 below. 


2. Boundary Conditions 
The PAL2 software accepts forces, pressure loads, line loads, and 
concentrated moments. This software limitation requires the boundary 
condition moments and shears to be distributed by means of a hand 
calculation. The calculated shear and moment loads apply indirectly to the 
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nodes on the boundary of interest via a dummy structure that when coupled 
with the Saint-Venant principle generates the centroidal load conditions on 
the center section boundaries. This last procedure is necessitated by the 
need to generate centroidal loads to validate the hand calculations. Once 
the model has been validated, non-centroidal loads may be applied to 
develop a complete structural response envelope given a variety of inertial 
load conditions. Once an inertial load is specified for a given flight 
condition, the structure is equilibrated; Displacements and rotations 
specified as zero at one node will prevent it from flying off into space. 


3. Model One 

The first model built contains the rod <nd shec: structure used in 
the hand calculation (see Figure 12). The model sizing comes directly 
from the hand calculation with one exception. The wing box thickness was 
modified to accommodate the lattice structure. A compromise height of 
sixteen inches required a slight adjustment (less than 1 part in 40) to the 
rod cross sectional areas. This adjustment maintained constant moment of 
inertia terms I, and Ij. The node locations/model dimensions are listed 


in Appendix C. The bar diameters were taken directly from the hand 
calculation Aj and A; terms except where they are altered to account for 


the change in wing box thickness. For simplicity, the shear plate thickness 
was uniformly chosen as 0.25 inches. 
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Figure 12. Center Section Wing Box Structure : model one (bars 
and sheets) 


Modeling the hand calculated centroidal load requires a dummy 
structure built from the face of the center section on each of its four 
boundaries. This dummy structure enables the centroidal loads to be 
applied at the boundary by means of the Saint-Venant principle in which 
only the nonequilibrated loads are transmitted through the structure. These 
nonequilibrated loads transform to the centroidal loads as they propagate 
away from the boundaries. 

At wing station 65 for example, the rectangular cross section consists 
of four nodes. This plane is repeated two times in both the positive and 
negative y-direction as depicted in Figure 13. There are a total of four 
sheets on the y faces and eight each on the x and z faces of the dummy 
structure. Rods connect each pair of adjacent nodes. The hand calculation 
sizes the rod diameter. The loads applied outboard of wing station 65 are 
such that at the boundary of interest they become the centroidal loads 
identical to those used in the hand calculation. The four nodes on the 
positive y face of the resulting structure have their displacements set to 
zero to ensure equilibrium throughout the structure. 
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WING STATION 65 


Figure 13. Wing Station 65 Boundary : modeled for comparison 
with hand calculations 


The remaining three boundaries are modele:. in a parallel manner 
using the dummy structure to achieve a centroidal loading on the boundary 
of interest for direct comparison to the hand results. 

The verifications of model formulation and implementations are made 
against a hand calculated bench mark model for the simplified centroidal 
loading conditions (described in Chapter II). The results are described in 
Appendix D. They are well within the theoretical expectations and the 
model presented herein can be considered as fully verified; it can be used 
for feasibility studies, operational envelope definition and structural design 
tradeoffs. 
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VI. ALTERNATIVE SOLUTIONS (PROBLEM EXTENSION) 


This chapter discusses the limitations of verifying the finite element 
solution with the hand calculated results to the boundaries of the structure 
for a 1g centroidally distributed load. Extensions from these verified 
boundary results explore the full range of the finite element solutions and 
their applicability to the existing center section wing box problem. The 
solution range encompasses all of the loads discussed in Chapter II. 


A. BENEFITS OF THE DUAL APPROACH 

The usefulness of the combined solution method falls into two 
categories. The first lies in the verification of the centroidal boundary 
condition (or lower bound) stresses for the finite element method. Ít 
should be pointed out that the results of the hand calculation are valid only 
for the boundaries of the model. This limitation lies in the generation of 
the boundary condition loads. The loads are generated external to the wing 
box and do not account for any distribution of inertia loads within the 
center section itself. Secondly, completion of the finite element modeling 
enables a rapid solution to the altered boundary conditions associated with a 
3g load for a 95,000 pound zero fuel weight P-3 with 60,000 pounds of 
fuel (the original configuration of interest). Since the solution varies 
within the center section as the inertial loads are redistributed, several may 
be applied simulating any of a number of configuration options. The most 
straightforward is an equitable distribution throughout the center section. 
From this original problem, design options are explored which will satisfy 
the structural requirements associated with the proposed payload growth. 


1. Configuration Change 
A configuration change to the P-3H short of a structural 
modification requires that a solution to the altered boundary condition 
loads be applied to the existing finite element model. The internal stresses 
are solved using the available software. The new fuselage boundary 
condition loads are calculated in the same manner as those calculated in 
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Appendix A. For the case of incre: ing the zero fuel weight, additional lift 
must be generated by the wings to „stain this extra mass at the 3g flight 
condition. The second solution found in Chapter II under wing boundary 
condition loads, provides a component build up method to calculate the new 
moments and shears at wing station 65. 


a. Example 1 

suppose a new ESM suite at SS3 adds 100 pounds to the work 
station 100 inches forward of F.S. 571. The 3g flight condition requires an 
additional 300 pounds of lift out of the wings distributed along the span. 
This addition alters the moment and shears at wing stat' 1 65. Fuselage 
station 571 would experience an additional 300 pounds of -hear and 30,000 
in-pounds of bending moment. These new loads are added to the existing 
boundary condition loads and a solution generated. The internal stresses 
may then be examined for failure determination. In the case where no 
overstress has occurred then the inertia loads within the center section may 
be redistributed until the limiting stresses are reached. 


2. Changing Material Properties 

There exist three solutions to the problem load which exceeds the 
yield stress of a given structure given constant center section inertia forces. 
The first two alternatives increase the moment of inertia terms 
(denominator of Equations 3.21-3.26) the `y reducing the internal stress. 
Either additional material may be used к acrease the cross section (ie. 
fortify existing sections) or the dimensions of the wing box and/or 
fuselage must be increased. The more appealing third option from the 
weight viewpoint involves the use of a material with a greater specific 
Strength. 

Replacing the entire aluminum structure with a single material (i.e. 
composite) fails to alter the internal stresses. Rather, the yield stress of the 
structure increases enabling greater loads to be safely carried. The partial 
use of stronger materials in high stress locations of the wing box would 


solve the structural problem for the case in which only portions of the 
structure exceeded the yield stress. 
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a. Example 2 

Changing the load requirements from 3.0 to 3.5 g's provides another 
opportunity to observe the structural limits. The problem requires that the 
3.5 g boundary load conditions for the wing and fuselage stations be 
calculated. As a rudimentary approach, the 3 g loads are multiplied by the 
ratio (3.5/2.0). The new loads are applied to the boundaries of the 
computer model and the stresses checked. Exceeding the yield stress 
anywhere in the structure requires a redesign by one of the three methods 
mentioned above. Assuming that only the stress in the axial bars between 
wing station 65 and the fuselage are excessive they will be targeted with 
composite replacement. The aluminum will be replaced by a composite 
material with a simple software modification. The new model will be run 
against the 3.5g boundary condition loads to complete the process. The 
new internal stresses within the member can now be analyzed, and further 
modifications made if necessary. 

Two extensions from the original verified finite element solution are 
explored. The first discusses the problems associated with payload increase 
consisting of a single piece of equipment, while the latter confronts the 
issues surrounding an increase in the flight envelope of the aircraft. Both 
extensions demonstrate design benefits of the verified finite element 
solution for the purpose of answering the structural integrity question 
associated with the payload increase. 
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УП. RESULTS 


This chapter examines the internal stresses of the P-3C center section 
wing box for the two calculation methods that were employed in the 
analysis. The hand calculations are based upon the Euler-Bernoulli 
assumption that plane sections remain plane. The finite element method 
simultaneously satisfies the field equations for a given element from which 
the internal stresses are determined. The results compare a 1g centroidal 
boundary load condition (ie. lower bound). 

Comparison of the stresses is limited to wing st^^ ^n 65L and fuselage 
station 571. Similar results are expected for the er two bounc es. 
The right wing boundary condition is symmetric to the left. Fuselage 
Station 695 has a scaled down shear force while its cross-section is identical 
to the forward section. . 

The results are presented graphically (Figures 16 and 17) and followed 
in tabular form (Tables 7.1 through 7.4). Figures 14 and 15 depict the 
internal stress locations for F.S.571 (X-face) and W.S.65L (Y -face). 


Ж 
Wing Station 


65 


Fuselage Station 
571 







Figure 14. Stress Locations 
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Figure 15. Stress Identification : X face(top) Y face(bottom) 
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Figure 16. Graph of Shear Stress Results 
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Figure 17. Graph of Axial Stress Results 
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TABLE 7.1. X Face Axial Stress 


FEM (psi) Hand Calc (psi) 
1159 1176 
-5.5 11.8 
-294 -278 
-333 -314 
-4.9 11.8 
-294 -278 


TAB: 7.2. X Face Shear Stress 


FEM (psi) Hand Calc (psi) 

0 12 
313 315 
0 24 

0 16 

0 15 
108 95 
| 165 
i 100 
12-5 21 
28 25 


TABLE 7.3. Y Face Axial Stress 


FEM (psi) Hand Calc (ps1) 
8932 9355 
-8932 -9355 
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TABLE 7.4. Y Face Shear Stress 


Shear Stress FEM (psi) Hand Calc (psi) 
S 3361 3568 
T -569 -518.5 
U 569 518.5 
V 2369 2551 


The stresses on the X-face of the model result from an application of a 
12,566 pound shear force 260 inches from the boundary of interest. This 
load results in the correct shear and moment at F.S.571 due to the Saint- 
Venant effect. With this centroidal loading of the fuselage, the stresses in 
the horizontal sheets are expected to be zero. The hand calculation leaves 
small finite stresses due mainly to round off in the five by five matrix 
solution. The finite element method predicts the stress as zero. The finite 
element axial stress in member A tends toward the high side, while the 
remaining rod stresses tend toward the low. This phenomenon results 
from the centroidal axis of the finite element model shifting towards the 
top of the fuselage since the area of the shear sheets are included. 

The shear load applied at thirty-three percent chord combined with the 
force couple realistically models the aerodynamic loads on the wing. The 
finite element stresses on this face are all within ten percent of the hand 
calculated results. 

The finite element stresses agree favorably with the hand results, well 
within twenty percent. Axial stresses tended to be closer to the predicted 
hand results than did the shear stresses. The Euler-Bernoulli assumption 
explains the difference between the two sets of results. 

These results verify the finite element model. The extension from this 
problem remains to run the model using a realistic 3g load condition for a 
95,000 pound zero fuel weight aircraft with 60,000 pounds of fuel. Since 
the fuel and internal loadout are never constant, the inertia loads within the 
center section must be varied to develop a full range of solutions for the 
center section wing box. 
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VIIT. SUMMARY, CONCLUSIONS AND 
RECOMMENDATIONS 


The Navy 1s considering the feasibility of increasing the patrol aircraft 
P-3C zero fuel weight enabling avionics and payload growth. This analysis 
examines the consequences to the structural requirements of the center 
section wing box. Specifically, the investigation provides an analytical 
means with which to explore the structure given a payload increase of 
18,000 pounds. 

Two solutions to the structures field equations are investigated: a 
simplified hand solution for preliminary feasibility calculations and a more 
precise solution for design analysis. Together, the solutions provide a 
necessary check for the results. The hand calculations calibrate the finite 
element solution on the boundary of interest. These calculations are valid 
only for the extremely simplified centroidal loading. This loading gives 
rise to uniformly distributed forces and moments and hence, the best case 
(lower bound stress) condition. The primary purpose of the hand 
calculation is to identify the existence of any structural problems on the 
wing box boundary and to assure the correctness of the finite element 
method model. This assurance is especially important wher the computer 
generated model is used to explore different options of payload 
redistribution feasibility studies, operation envelope definition and 
structural design tradeoffs. 

The extension remains to utilize the finite element model for various 
load conditions combining the boundary shears and moments together with 
several inertial distributions representative of the limiting flight envelope. 
The critical center section components are to be identified and redesigned 
from the wing box feasibility studies which, in turn, will lead to an 
optimized solution. The options will include structural fortification, 
composite reinforcement, substitution, and inertial distribution limitations 
within the center section. 
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At this stage, finite element output is limited to stress and deformation 
data. Software options are also available to generate dynamic responses 
which are needed for flutter characterization. 

This analytical package facilitates the generation of design points for 
the replacement center section wing box. Successful implementation of the 
product will provide the technical basis for RFP generation. 
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APPENDIX A. BOUNDARY CONDITIONS 


A. P-3C FUSELAGE WEIGHT BREAKDOWN 
(Weights from P-3C specification 1983) 


COMPONENT 


TAIL 
WING 


PROPULSION 
NACELLE 
LANDING GEAR 


AIR COND./ANTT-ICE 


FLIGHT CONTROLS 


ELECTRICAL 


INSTRUMENTS 
HYDRAULIC 
BODY 


WING 
C.S. 


MAIN 
NOSE 


WING 
NOSE 
C.S. 


WING 
C.S. 


WING 
CS: 


NOSE 
C.S. 
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WEIGHT(LBS) 
1965 
9095 
D95 
1500 
15041 
4877 
3666 
3055 
611 
2022 
800 
522 
800 
100 
1509 
750 
600 
Ey 
1813 
500 
1313 
477 
433 
9973 
4773 
1200 
4000 


APU 
ELECTRONICS GROUP 


ARMAMENT 
FURNISHINGS 


PHOTO 


NOSE 
( 


NOSE 
CS 


514 


10376 


1404 
3830 


89 
67,084 


B. SHEAR AND MOMENTS BY SECTION 


1. Nose -- 


0000000000000 


СОМРОМЕМТ 


INSTRUMENTS 
BODY 

APU 

ELECTRONICS GP 
ARMAMENT 
FURNISHINGS 

AIR COND. /ANTI-ICE 
PHOTO 


NOSE GEAR 
SHEAR = 12566 LBS 
SHEAR = 37698 LBS 


Loads at 


(LBS) 


477 
4773 
514 
33716 
404 
2000 
322 
89 
611 


MOMENT = 3.31 * 106 IN-LBS 


F.S. 571 


350 
275 
283 
200 
100 
300 
470 
470 
300 


3376 
6000 
1000 


2000 
330 
1500 


WEIGHT MOMENT-ARM MOMENT 
(IN) 


(105 IN-LB) 


1.66 
13.10 
].45 
6.75 
0.40 
6.00 
ЕЭ 
0.42 
1.83 


@ lg 


MOMENT = 9.93 * 10© IN-LBS @ 3g 
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2. Rear -- Loads at F.S. 695 


COMPONENT  ............. WEIGHT MOMENT-ARM MOMENT 
(LBS) (IN) (10? IN-LB) 
TAIL GP 1965 505 9.92 
BODY 4000 250 10.00 
FURNISHING 1500 200 3.00 
FLIGHT CONTROLS 159 300 0.48 
ELECTRONICS 1000 100 1.00 
AIR COND. /ANTI-ICE 100 400 0.40 
SHEA: = 8724 LBS {OMENT = 2.48 * 106 IN-LBs @ Ig 
SHEAR = 26172 LBS MOMENT = 7.44 * 106 IN-LBS @ 3р 


EQUILIBRIUM REQUIRES NOSE & AFT MOMENT EQUAL 
THEREFORE ASSUME DIFFERENCE IN MOMENTS COMES FROM 
ELEVATOR FORCE. 


ELEVATOR MOMENT = 0.83 * 106 IN-LBS 


3. C -nter Section 


COMPONENT WEIGHT (L.5S) 
HYDRAULIC 433 
ELECTRICAL 1515 
FLIGHT CONTROLS 600 
BODY 1200 
ELECTRONICS GP 6000 
ARMAMENT 1000 
FURNISHINGS 330 
AIR COND. /ANTI-ICE 800 
WING 1500 


CENTER SECTION WEIGHT = 13176 LBS @ lg, 39528 LBS @ 3g 
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C. WING LOADS 
(Provided by Aerostructures Inc.) 


200000 


150000 


100000 


Shear (Ibs) 


50000 





Wing Station (in) 


Figure 18. P-3C Wing Shear : 135,000 pounds gross weight. 
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Moment (in-Ibs E404) 


Lnd | IA E P 





Wing Station (in) 


Figure 19. P-3C Wing Moment : 
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135,000 pounds gross weight. 


APPENDIX B. GENERAL SHEAR SOLUTION 


The following argument follows directly from Reference 1. 


` q 
— == б 
XS 
| t 
X 
IG ) 
«— 4— <q <q 


— — — — 
q(x,s1) +A QE 


с o‏ )ل هه 


Solving for the shear begins with summing the forces in the x 
direction: 


YF,-0 (BĲ) 


s (r^s s r^s 
0 - | O,, t(Xos) ds - | (GAG, [t(xos)- At,(s)] ds + 


1 1 


J 4(х,51) +Aq,) dx - J q(x, sS) dx (B.2) 


xO-Ax 
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Cancelling like terms and dividing by AsAx Equation B.2 simplifies to: 





1 At iE 

= б, go ds + 

AS 

xÜ 

1 Aq (x) 1 At (s) 
— — dx- — Ao... — ds (B.3) 
Ах AS AS * Ax 

x0-Ax s 


Since the thickness does not vary along the length į: the sheet those te: 
fall out of the expression leaving: 


дс 
е -t(x ps) — 
ds 


(B.4) 





Integrating both sides of of Equation B.4 yields: 


S1 


Oxx 


9 
q(s,) = q(s=0) - t = ds (В.5) 





X 


0 


Substituting the general expression for axial stress (6,,) developed іп the 


earlier solution into Equation B.5 produces: 
І 


M M 

q(s;) = q(s=0) - Si Ba Mt Ml, Mon Eten, tds (B.6) 

Әх E. * E, * * 
Cyan” í) (yy 7 Lye) 


0 


Applying symmetr to the geometry removes the cross product of inertia 
terms. Taking the: tial derivative with respect to x yields: 


D 


S1 





М.Е І] ðM, E 1 
q(s ) = q(s=0) - { —— H ` Ve 1 | J: ds (B.7) 
dx E, I dx E, I 
ZZ yy 
ӘМ, 


OM 
Dor Ve (6% "s -m (x) +V (x) (B.9) 
X X 








Substituting the definitions from Equation B.8 &9 into B.7 and neglecting 
the distributed moments for the given 3g loading condition results in: 





Э] 
Е -V,@) Е Ух) 
q(s;) - 4(5-0) - { —{——Jy +4 : Jz}tds (B.10) 
Pa Iz š lyy 


For a given location of x the shears are constant and may be removed from 
the integral along with the moments of inertia. Employing the definitions 
from Equation B.11&12, Equation B.10 reduces to: 


51 51 
ж 


Е ў E 
Q= Í y tds (B11), Q= f ата (B.12) 
E, E, 


V (x) уху» 
q(s;)= q(s=0) - ——Q,-——Q, (B.13) 


p 1 vy 


Assuming the shear (Vy) produced by a rudder deflection and experienced 


at F.S. 695 1s negligible, equation B.13 reduces to: 


Vx) * 
q(sp-q(s-0 -——Q, (В.14) 


lyy 
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APPENDIX C. NODAL POINT LOCATIONS FOR FEM 


A. FEM OF FUSELAGE BOUNDARY CONDITION 


П 


Fuselage Boundary Modeling 
nodal point locations 1 

1,0, 54, -72 through 5, 527, 54, -72 
6, 0, -54, -72 through 10, 527, -54, -72 
11,0, 54, -8 through 15, 527, 54, -8 
16, 0, -54, -8 through 20, 527, -54, -8 
21,0, 54, 8 through 25, 527, 54, 8 
26, 0, -54, 8 through 30, 527, -54, 8 
31,0, 54, 10 through 35, 527, 54, 10 
36, 0, -54, 10 through 40, 527, -54, 10 
41, 0, -63, -8 through 45, 527, -63, -8 
46, 0, -63, 8 through 50, 527, -63, 8 
51,0, 63, -8 through 55, 527, 63, -8 
56,0, 63, 8 through 60, 527, 63, 8 
blank line 
C bars are all circular cross-sections 
c quad plates are all shear members 1/4 inch thick 
c 2024-T4 aluminum material 
material properties 10.5E6, 0, 2.59E-4, 0.33,40E3 
c X dir rods 
beam type 3 4.76 0 
element generate 21 

1510150 
blank line 
Беат type 3 6.31 0 
element generate 21 
5155 45 1100 
56 60 501 100 
212530150 
blank line 
beam type 3 7.90 0 
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element generate 21 
[1 1520150 

blank line 

beam type 3 2.52 0 
element generate 21 
31 35 40 1 50 

blank line 

c Y dir rods 

beam type 3 2.12 0 
element generate 21 
1510151 
313540151 

blank line 

beamtype 3 5.15 O 
element generate 21 
45 20 162510 
201511510 
1555514010 
5030262010 
302521510 

25 60 56 35 10 
blank line 

c  Z dir rods 

beam type3 2.12 O0 
element generate 31 

1 510401 5 102 
blank line 

beam type3 1.78 0 
element generate 31 
41 45 55 601 1052 
blank line 

quad plate element 3 0 0.25 
element generate 31 
15104015100 
15104015101 
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1 5 10 40 1 5 10 2 

41 45 55 60 1 10 5 1 

blank line 

do connect 26 46 41 16 through 30 50 45 20 step 111 1 
do connect 21 56 51 11 through 25 60 55 15 step 1 1 1 1 
element generate 21 

16 20 45 1 25 

26 30 50 1 20 

5] 55 15 1 40 

56 60 25 1 35 

blank line 

end definition 


2. Fuselage Load Conditions 
c load condition for fuselage station 571-lg St Venant 
forces and moments applied 1 
fz 2094 51015 20 25 30 
blank line 
displacements applied 1 
ta 0.0 16 
tx 0.0 1611 21 2631 36 41 4651 56 
tz 0.0 1611 21 2631 36 41 4651 56 
blank line 
solve 
quit 


B. FEM OF WING BOUNDARY CONDITIONS 
1. Wing Boundary Modeling 

c wing station 65 shear loads located at 33% chord 
nodal point locations 1 

1 0-60 0 

2 80-60 0 

3 0-60 16 

4 80 -60 16 
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17 0-70 0 

18 80-70 0 

19 0-7016 

20 80-70 16 

blank line 

nodal point locations 31 

12420124 

blank line 

c all beams have circular cross section 
c quad plates are all shear members 1/4 inch thick 
c 2024-T4 aluminum material 
material properties 10.5E+06, 0, 2.59E-4, 0.33, 40E+03 
beam type 3 5.04 0 

element generate 31 

1242012420 

blank line 

beam type 3 1.38 0 

element generate 31 

1242012400 

1242012410 

blank line 

quad plate element 3 0 0.25 

element generate 31 

124201240 

124201241 

124201242 

blank line 

end definition 


Wing Load Conditions 
c wing box centroidal loading lift line at 33% chord 
forces and moments applied 1 
fy 1.828E+05 19 20 
fy -1.828E+05 1718 
fz -5064 17 19 
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fz 7136 18 20 
fx 5184 17 18 
fx -5184 19 20 
blank line 

displacements applied 1 

ta 0.01234 

blank line 

solve 

quit 
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APPENDIX D. FINITE ELEMENT IMPLEMENTATIONS 


Appendix D outlines the procedures necessary to correctly model a 
structure using the finite element method. The first portion discusses the 
modeling elements and is followed by an explanation of displacement and 
equilibrium force boundary conditions. Application of the Saint-Venant 
Principle to modeling centroidal force distributions completes the topic. 
Examples which highlight these issues are supported with numerical results 
and summarized. 

The first decision that one must make when modeling a structure lies in 
choosing the correct elements. Rods and sheets constitute the model for the 
hand calculation. The PAL2 software provides four beam type elements 
and two quadrilateral plate elements from which to select. The first sheet 
element provides resistance to bending, axial and shear loads. In resisting 
bending and axial loads, this element impedes relative displacement motion 
between the nodes that surround it. Another way to state this behavior 1s 
that the element stiffens the k matrix. The second sheet type resists only 
shear loads. Since the plate elements do not resist relative motion between 
the surrounding nodes, this plate must be surrounded by beam elements as 
a picture is surrounded on its four sides by a frame. Beam 4 models an 
anti-symmetric cross sectional element. Beam 3 models a rod with circular 
cross section. The user defines the inner and outer rod diameter. The 
program automatically calculates the remaining geometric properties from 
this information including the shear area. The shear area represents a 
fictitious area that corrects for the parabolic shear distribution across the 
face of an element. For a circular rod, the shear area equals the cross 
sectional area divided by 1.185. Beam 2 creates a rectangular cross section 
which varies in height and width along its length. Beam 1 generates a 
constant area cross section member. In addition to the area, the user must 
specify the torsional moment of inertia, the moments of inertia, the shear 
area, and the distance from centroid to the most distant point on the cross 
sectional area. The shear area correction factor for a square cross section 
beam is 1.20. Curved beam elements are also available but not used since 
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the hand model was constructed from straight rods to facilitate calculation. 
The second area of interest lies in displacement boundary conditions; 
how to tie the structure down in space to prevent that unequilibrated shear 


of 1077 pounds from driving the structure into the computer abyss. If six 
degrees of freedom were eliminated from one location on the structure, 
would the tie down affect the internal stresses? To illustrate this point, 
four circular rods of one square inch cross section (beam type 3) surround 
a 100 square inch shear plate of the second type discussed above. The left 
two nodes have their x, y and z translations set to zero while a shear force 
of unity in the negative z direction is applied to the upper right hand node. 
The reaction forces are calculated at the left end by the software. “ese 
reaction forces are applied to the left end with the same shear load . hile 
the displacement boundary conditions are eliminated from the left end and 
three translations and rotations are set to zero on the lower right hand node 
of the structure. The reaction forces calculated at the lower right hand 
node for this second case are eleven orders of magnitude less than the 
applied load while the internal stresses are unchanged in all five elements. 
similar results are obtained for an example extended into the third 
dimension. The results from this exercise enable the center section wing 
box to be tied down at one node since it too is in static equilibrium. With 
confidence one knows that there will be insignificant stresses generated 
within the modeled structure due to the displacement boundary condition. 
Finally, the St. Venar principle must : addressed when applying 
boundary condition loads and displacements. Unequilibrated loads are 
transmitted through the body while equilibrated loads dissipate quickly as 
they extend through the structure. From the modeling aspect this affects 
the manner in which one exerts forces on the structure. If large single 
point forces are applied at a few nodes around the structure, this practice 
may lead to the finite element solution being driven to an alternate form. 
Consider a hollow square beam forty inches in length (figure 20). 
The beam consists of four cubes of eight nodes each whose faces in the 
longitudinal direction are absent while the remaining four faces are made 
from shear plates (type 2) 1/10 inch thick. Circular rods of five square 
inch cross sections connect the longitudinal edges of the forty inc! eam. 
Square beams of 1/100 square inch cross sections surround the rez. ining 
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sides of the quadrilateral plate elements. The small dimension of the 
secondary beams enable a blending with the sheets. The four left end nodes 
of the long beam are attached with three translational degrees of freedom 
set to zero. The right end has a shear load of unity applied to the upper 
node nearest the observer. The manner in which one applies the load 
determines the form of the solution taken by the finite element method. If 
the load is applied as described above, the right end panel nearest the 
observer will display a shear of 0.984 psi. Observation of the output 
stresses five inches from the right end of the long beam reveals that the 
beam 1s in static equilibrium. Results from the hand calculation predict a 
value of 0.750 psi.(see Allen and Haisler p221). However, if the shear load 
above is redistributed to the centroid as a shear and moment with one 
quarter of each divided among the four end nodes, the expected solution of 
0.742 psi results or one percent variation from the hand calculated 
solution. Static equilibrium is satisfied by this solution at the same location 
described above. In practice, loads should be distributed as evenly about 
the structure as possible. 
2 





Figure 20. Square Beam Model 
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For the case in which the boundary condition loads are not perfectly 
matched to the internal stresses, Saint-Venant points out that the structure 
requires space to enable the equilibrated forces to dissipate. As an 
extension of the previous problem, consider a similar beam consisting of 
twelve cubes of which five are shown in Figure 21. The first plane of 
interest lies 1/8 th of the distance from the loaded end to the wall. The end 
cube of the beam is subdivided from one to two sheets per face (case a), 
from two to four (case b), and from four to six (case c). The increased 
number of end nodes facilitates the application of a more evenly distributed 
superposition of shear and moment loads to the beam in comparison to the 
previous example. The results of the three plar are compared :^ the 
hand calculations for each of the three conditions shown in Figure 21. 
Figure 22 shows the loading for cases a and b. Case c distributes the load 
in a similar fashion among sixteen nodes. 


Z 
X 
case a 


x case b 


Figure 21. Twelve Cube Beam 
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moment = 1/12 
shear = 1/8 + “p.p. 


+ v + 
I + 

Y 

I ШЕ 
t Y 


4 4 4 
Figure 22. Load Conditions : superposition of shear and torsion loads 
as seen from end view of case a and case b. 


A. SUMMARY 
In Figure 23 the observer is viewed from the negative z axis from the 
wall. In Figures 24, 25, and 26, Tj is located in the upper right-hand 


corner, -T9 = T4 is in the lower left corner, and T3 is in the center. 


The stresses in plane 1 (Figure 24) vary by as much as forty percent 
from the hand results while by plane 2 (Figure.25) the variation drops 
dramatically to 5.6 percent. The stress difference in plane 3 (Figure.26) 
reduces further to 3.2 percent. As the distance from the boundary 
condition increases from plane 1 to plane 3, the finite element stresses 
approach the hand calculations thereby validating the Saint-Venant 
principle. 


63 


PLANE 1 DATA 


2 
е p. 
PLANE OF Y 
T 
INTEREST Š 
1; | 


Figure 23. Axial View of Example Beams 


-0.05 0.15 0.35 0.55 
HAND 


CALCULATIONS 


Figure 24. Plane 1 Results 
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PLANE 2 DATA 


PLANE 3 DATA 


-0.25 





HAND 
CALCULATIONS 


Figure.25. Plane 2 Results 


ШІЛ 
к 
TEE 

Cer 


-0.05 


-0.25 -0.05 0.15 HAND 0.35 0.55 0.75 
CALCULATIONS 


Figure 26. Plane 3 Results 
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As an extension of these example problems to the center section wing 
box, dummy structural area must be located between the boundaries of the 
structure and the applied loading conditions. This area will ensure that the 
proper shear stresses are retrieved from the structure as occurred in 
moving from plane 1 to plane 2 above. The moments of the applied 
loading condition must be modified to account for the added distance 
through which they must act. 

These results which compare hand calculations (using the Euler- 
Bernoulli assumption) to those of the finite element method illustrate the 
importance of correctly modeling the boundary conditions to minimize 
their effect on the results. In addition, "dummy extended space" facilitates 
the location of centroidally distributed loads through the Saint-Venant 
Principle. 
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